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We discuss some applications of higher symmetry groups to condensed matter 
systems. We give special attention to the groups SO{n) (n = 4,5,6,8) in the 
two-dimensional Hubbard model and its generalizations, which model the high Tc 
cuprate superconductors. 



This talk is intended to be a brief introduction to the interesting roles 
played by higher symmetries in the theory of the cuprate high Tc supercon- 
ductors and other quasi-two-dimensional condensed matter systems. Since this 
in mainly a particle physics audience, I will give a quick introduction to the 
systems we are dealing with. A broad modern survey of superconductivity is 
given by Tinkhamlll, for example. 

1 High Tc superconductors 

The typical cuprate superconductor has Cu02 layers sandwiched between var- 
ious layers, which serve as charge reservoirs to provide carriers (holes, in fact) 
to the conducting Cu02 planes, as shown in fig. |l|. A typical charge reservoir 
is La2~xSrx02, in which replacing La^^ by Sr^ leads to a hole density x 
in the conducting layers. High Tc is typically obtained for doping fractions 
X ^ 0.15 — 0.2. Since superconductivity is mainly in the Cu02 planes, it is 
useful to treat the superconducting system as two-dimensional, with effects 
due to coupling between the layers ignored at first. 

2 Two-dimensional lattice 

A realistic model Hamiltonian for the Cu02 planes involves electrons (or holes, 
as is the case in practice) hopping between p-shell orbitals in O and c?-shell 
orbitals in Cu together with interaction terms which model the Coulomb repul- 
sion between two electrons on the same atom, or on adjacent Cil and O atoms. 
Detailed models ar£ reviewed by Kresin, Morawitz and WolfH. However, it 
has been suggestedcl that much of the essential physics of the Cu02 planes is 
captured by a simpler model in which holes only hop between Cu sites; a.more 
detailed exposition of this view is given in the recent work of Andersoncl. 
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Figure 1: Schematic views of the cuprate high Tc superconductors. The left figure (a) shows 
the conducting Cu02 layers (white) between the cross-hatched layers which serve as charge 
reservoirs, as explained in the text. The right figure (b) shows a typical Cu02 layer, with 
Cu atoms (dark circles) at the vertices of a square lattice and O atoms (open circles) at the 
midpoints of the lines joining adjacent Cu atoms. 



Thus we consider first the liopping Hamiltonian 

Ho = -t aLflya 4aaya 

<xy>,Q «xy»,a 

Here < > (« ») denotes sum over (next) nearest neighbor pairs, and a is a 
spin index. The spectrum of Hq consists of a single energy band, with energies 

ii^k — — 2t(cos kx + cos ky) ~ At' cos k^ cos ky (2) 

as k runs over the first Brillouin zone, which is depicted in fig. |^(a). At the 
points indicated by large circles in the figure, we have Vk-E = 0, and hence 
a singularity in the density of states, the Van Hove singularity (VHS), which 
is logapithmic in two dimensions. One of us (RSM) has written an extensive 
reviewB describing the importance of the VHS in high physics. 

Another feature important in understanding this lattice is the existence of 
a nesting vector Q such that exp(iQ • x) = ±1 for every lattice site x. This 
vector connects the Van Hove singularities, and it also joins points near the 
Fermi surface for electron (or hole) densities near one per site. This creates 
the possibility of various instabilities associated with wave vector Q, including 
antiferromagnetism, d-wave superconductivity, and others. Here the vector 
Q = (±7r, ±7r) when length is expressed in units of the lattice spacing. 

TcLdescribe the interactions of the electrons, we can start with the Hubbard 
modelB interaction 

Vnubbard = 51 "^xT^xi (3) 

X 

which represents the Coulomb interaction between two electrons on the same 
site. SchulzB has considered additional interaction terms. 
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Figure 2: First Brillouin zone for the square lattice, (a) The dashed hues show the Fermi 
surface for the simplest nearest neighbor hopping Hamiltonian (hi) with t' = and one 
electron per site {half filling of the band). The large circles at k = (0, itTr), (±7r, 0) give the 
location of the Van Hove singularities in the density of states (see text), (b) Division into 
two subzones X and y-. X contains the unshaded region together with the line kx + ky = 
and y contains the shaded region and the dark line kx = ky. 



3 Pairing algebras 

Consider fermion creation and annihilation operators a\,ai {i — 1, . . . , N). The 
N'^ number conserving operators alaj generate an SU{N) U{1) Lie algebra, 
and when we combine them with the N(N — 1) pairing operators 

Ujk = a.ftj and n]^^ ee a]a| (4) 

we obtain a pairing algebra which generates a pairing group SO{2N). Some 
elements of this algebra may be conserved; other elements O may satisfy a 
commutation relation of the form 

[H, O] = AO (5) 

with the Hamiltonian H. These elements define a spectrum generating algebra, 
since for an operator O which satisfies (H), 

Hi;^ EtP^ H{Oip) = {E + X) (6) 

Thus the operator O transforms one energy eigenstate into another with eigen- 
value different by A. The first systematic study of these algebras iujcondensed 
matter physics appears to have been done by Birman and Solomonla. 
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4 A simplified SO (8) algebra 

A simplified 5*0(8) algebra can be constructed if we let aj^,, be creation and 
annihilation operators for a fermion with momentum near (±7r, 0) (a, (3 are 
now spin indices), and the corresponding operators for momentum near 

(0, ±7r). This choice is motivated by the existence of the Van Hove singularities 
in the density of states, which lead us to suspect that electrons (or holes) 
with momenta near these singularities will be most important (since there are 
more of them). As explained in the previous section, we obtain from these 
operators an SU{A) © C/(l) algebra of number conserving operators and an 
50(8) pairing algebra, some elements of which were used by SchulaJ to analyze 
pairing instabilities in generalized Hubbard models. 
The operators 

S = ^ {ai~^apap + bi~^^pbp^ ' ^ " \ ("""^"/5"/3 ~ (7) 

and 

{aW + bib„ - 1) (8) 

generate an 50(4) © C/(l) algebra which is a symmetry of the nearest neighbor 
hopping Hamiltonian (yj) (with t' = 0). The remaining number conserving 
operators 

OcDW = l{aib^ + bia^) , Ojc = {aib^ - bU^) (10) 

(11) 

generate instabilities (charge density wave, spin density wave, etc.). 
The pairing operators 

77 = Qaba , H = aa~a^pbp (12) 

[(ia = ap{ia''^)pa] belong to the 50(6) spectrum generating algebra of the 
nearest neighbor Hamiltonian, and the additional pairing operators 

As=aaaa+baba , Ad = aatta - baba (13) 

in the 50(8) pairing algebra generate instabilities corresponding to s-wave and 
d-wave superconductivity, respectively. 
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5 The full SO (8) algebra 

We now describe the full 50(8) algebra constructed using the nesting vector 
Q introduced above. We start with the usual charge and spin operators 



Q^Y. - ^ ) = E (4„Cka - I) (14) 

^ = ^ E ''i" ^a/jflx/S ^IY. ^k" ^a/3Ck/3 (15) 

(we use o^, a to denote coordinate space operators, c^, c to denote momentum 
space operators). Now divide the Brillouin zone for the square lattice into 
two subzones X and y, as shown in fig. ||(b), and introduce the characteristic 
function (also used by HenleyS) 

_r+i, hex 
^^-{-i, key (1^) 

The operators T='=, T^ defined by 

- E E CkaCk+Q,a ' T- = ^ ^ 4aCk+Q,a (1^) 
keX a key a 

- ^ E -k4.Ck„ (18) 

k.a 

generate another SU{2) ~ 5*0(3) algebra, which we call isospin; the 'flavor' 
here is simply the direction of the largest momentum component. Also define 



Ti = i (T , . , 

2 ^ ^ 2 ^ 



T-) = ^ E ^o^k+Q,. (19) 
4 (T+ - T-) = E -kcLck+Q.. (20) 



2i ^ ^ 2i 

k,a 



The spin SU{2) and 'isospin' SU{2) can be enlarged to SU{A) in the manner 
of the supermultiplet models in nuclear theory by including the operators 

T+ = E E 4a^a;3Ck+Q,/3 ^ T- = ^ 5] cUc.pC^+QJ3 (21) 
keX Q/3 key a/3 

= ^ E ^kcL^a/3Ck/3 (22) 
k,Q/3 



and, analogous to ( p^ ) and (pO|), 

T^ = i(^T+ + T") , t2 = ^(^T+-T"^ (23) 

The number conserving SU{2) © SU{2) ~ 5*0(4) algebra generated by the 
operators S and T ^ commutes with the free particle Hamiltonian Hq of eq. (Q) 
[the commuting SU{2) generators are actually S j- = S ± T^, which are the 
spin generators restricted to the subzones -^(4-) and 3^(— )]• The operators Q, 

also commute with Hq. 

We can also introduce pairing operators 

A = - ^ CkaCkQ , At = 2 X! '^kCkaCkQ (24) 

and extended pairing operators 

n=- ^Ck+Q,aCkQ , n = - ^ TkCk+Q,Q CrQ/3Ck/3 (25) 
k.a k,Q/3 

[here Cka = (*CT^)Q/3Ck,/3 is a 'charge conjugate' operator which has the same 
transformation properties as the creation operator for an antiparticle of mo- 
mentum k, spin state a]. The A and 11 operators, together with their adjoints, 
augment the SU{4) ® J7(l) algebra of number-conserving operators to form an 
50(8) pairing algebra. The operators introduced here are equivalent to those 
introduced in the simplified version, as shown in table |l|. 

The largest subalgebra of the 5*0(8) which commutes with the simplest 
hopping Hamiltonian Hq [eq. (0) with t' = 0] is the 50(6) © 50(2) algebra 

generated by Q, S, T^, H, nt, H, nt and T^ [T^ is the 50(2) genera- 
tor]. Note that the subalgebra can also be described as SU{4:) © t/(l), since 
5'C/(4) ~ 5*0(6) and U{1) ~ 50(2) (this equivalence of Lie algebras is some- 
times overlooked in the literature). This 5*0(6) © 5*0(2) provides a starting 
point for analyzing collective instabilities. 



Table 1: Translation table for elements of SO(8) in the simplified form of sec. 4 and in the 
full form of sec. 5. The operators Q, ~S and Tf are common to both presentations. 
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6 Symmetries and instabilities 

The generators of 50(8) not in the SO{6) © 50(2) subalgebra can be grouped 
into two 6-dimensional superspins: 

I]± = (T±,T±,ATA,,At±A|) (26) 

which are eigenstates of the 50(2) generator T'^ with eigenvalues ±1. Gener- 
ahzed pairing operators of the form 

= ^ E f(k)CkaCkc. (27) 
k,a 

which generate BCS-hke states can be introduced. These can also be grouped 
into multiplets which transform like vectors under the 50(6) group, and de- 
fine instabilities which will compete when interactions which break the 50(6) 
symmetry paxe included. 

Zhang^il has considered the 50(5) algebra obtained from 50(6) by deleting 

the generators T^, 11, 11^ The 50(3) 'isospin' algebra generated by T±, T^ 
commutes with Zhang's 50(5), and we are left with a triplet of 5-dimensional 
superspins: 

Z± = (T±,ATA„At±At) , Z° = (T3,n,nt) (28) 

These algebras are not exact symmetries of any laaturar two-dimensional 
model (though an ad hoc model has been constructedi3) . The only known fully 
two dimensional model with a higher symmetry is the Hubbard model, with 
nearest neighbor hopping and interaction given by eo— /^). This model has an 
50(4) ~ SU(2) © SU (2) spectrum generating algebraic consisting of the usual 
spin SU{2) and a second SU{2) generated by H, II'I' and Q. There are other 
symmetries for two-site models of the type discussed in sec. ^, which will be 
discussed elsewhere. 

7 Conclusions 

We have described some higher symmetries and approximate symmetries which 
may be useful in classifying the phase structure of quasi-two-dimensional sys- 
tems such as the cuprate high Tc superconductors. Dynamical calculations ace 
required for detailed description - a first step in this direction, as been takenij, 
using methods developed originally by Balseiro and FalicovHa. 
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